In this paper we present the results of the 3D modeling of convection under cratons of Central Asia. Numerical model is based on the "vector potential - vorticity" variables and the method of fractional steps. Computational results that demonstrate the structure of convective fluxes